Multilevel Linear Sampling Method for Inverse Scattering Problems

A novel multilevel algorithm is presented for implementing the widely used linear sampling method in inverse obstacle scattering problems. The new method is shown to possess asymptotically optimal computational complexity. For an n×n sampling mesh in R2 or an n×n×n sampling mesh in R3, the proposed algorithm requires one to solve only O(nN−1) far-field equations for a RN problem (N=2,3), and this is in sharp contrast to the original version of the method which needs to solve nN far-field equations. Numerical experiments are presented to illustrate the promising feature of the algorithm in significantly reducing the computational cost of the linear sampling method.